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Abstract 

A one-dimensional gas of electrons interacting with long-range Coulomb 
forces (V(r) ~ 1/r) is investigated. The excitation spectrum consists of sep- 
arate collective charge and spin modes, with the charge excitation energies 
in agreement with RPA calculations. For arbitrarily weak Coulomb repul- 
sion density correlations at wavevector \kp decay extremely slowly and are 
best described as those of a one-dimensional Wigner crystal. Pinning of the 
Wigner crystal then leads to the nonlinear transport properties characteristic 
of CDW. The results allow a consistent interpretation of the plasmon and spin 
excitations observed in one-dimensional semiconductor structures, and sug- 
gest an interpretation of some of the observed features in terms of "spinons" . 

A possible explanation for nonlinear transport phenomena is given. 
71.45.-d, 72.15.Nj, 71.28.+d 



1 



Typeset using REVTpjX 



The properties of models of one-dimensional interacting electrons have been studied in 
great detail. Examples are the so-called "g-ology" model of fermions moving in a continuum 
]1].|2|], or the one-dimensional Hubbard model ||[1],|4]] . In these models, one usually assumes 
short-ranged (effective) electron-electron interactions. The so-called "Luttinger liquid" || 
behaviour in this type of models is characterized by separation between spin and charge 
degrees of freedom and by power-law correlation functions, with interaction-dependent ex- 
ponents. Short-range interactions are a reasonable assumption for applications e.g. to 
quasi-one-dimensional conductors, where screeening between adjacent chains leads to ef- 
fectively short range interactions within one chain ||. However, the situation can be quite 
different if an isolated system of electrons moving in one dimension is considered. There 
then is no interchain screening, and the true long-range character of the Coulomb forces 
(V(r) = e 2 /r) needs to be taken into account. This appears to be the case, e.g., in certain 
one-dimensional semiconductor structures, where the effects of one-dimensional Coulomb 
forces have indeed been observed 0. 

The purpose of the present paper is to investigate the effects of the long-range Coulomb 
interaction in a one-dimensional model, using the bosonization method This allows in 
particular a rather straightforward and asymptotically (for low energies and wavenumbers) 
exact description of excitation spectra and correlation functions. The main conclusion is 
that the long-range force, even if it is very weak, leads to a state characterized by quasi- 
long-range order much closer to a one-dimensional Wigner crystal || than to an electron 
liquid. The calculations presented here provide a rather simple microscopic description 
of the Wigner crystal, a problem that in higher dimesnions has been difficult to treat by 
many-body techniques. 

I will start by considering the particular case of one-dimensional electrons with a linear 
energy-momentum relation interacting with long-range Coulomb forces, described by the 
Hamiltonian 

^ = II^[(^-M4, s a M + (- fc -^)&i,A, s ] + ^II^(9)p g p- (? + ^ s • (i) 

k,s <? 
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Here a k s (b k s ) creates a right- (left-) moving electron with momentum k and spin projection 
s, Vp is the Fermi velocity. In the interaction term p q = p a (J + p bq is the Fourier component 
of the total particle density, and V(q) is the Fourier transform of the interaction potential. 
In strictly one dimension, a 1/r Coulomb interaction does not have a Fourier transform 
because of the divergence for r — > 0, however, in a system of finite transverse dimension 



d, the singularity is cut off a r « d ||. Using the approximate form V(r) = e 2 /-*/r 2 + d 2 , 
one has V(q) = 2e 2 K (qd). Finally, the backward scattering term H bs describes processes 
where particle go from the right- to the left-moving branch and vice versa. This involves a 
nonsingular interaction matrix element at q ~ 2k f, called g\. 

The linear energy-momentum relation makes the model ([!]) exactly solvable using stan- 
dard bosonization methods. Moreover, at least for weak Coulomb interactions, when states 
near the Fermi energy play the major role, linearizing the spectrum is not expected to change 
the physics drastically, and one therefore expects that the model ([!]) correctly represents the 
low-energy physics even for more realistic bandstructures. The model can be easily solved 
introducing the phase fields 

M*) = - 7 E -e~ ip > + (p) + v-(p)) , (2) 

where v — p,a, and p r (jp) (a r (p)) are the usual charge (spin) density operators for right- 
(r = +) and left-(r = — ) going fermions. The Hamiltonian then decomposes into commuting 
parts for the charge and spin degrees of freedom. The charge part takes the simple quadratic 
form 



p 2tt 



J dx [n\l + ~g 1 )Tl 2 p + {l-~g 1 ){d x <P P f 



+ I dxdx'V(x - x')d x (j) p d x i(j)p (3) 
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Here U p is the momentum density conjugate to 4> p , and g\ = gi/ftirvp)- The Hamiltonian 
(0) is quadratic in the bosonic fields and therefore can be diagonalized straightforwardly. 
The elementary excitations then are found to be charge oscillations (plasmons), with energy- 
momentum relation 



= v F \q\ [(1 + sfO(l - ft + 2K(g))] 1/2 (4) 

where V^(g) = V A (g)/(7Tfir). The long-wavelength form, 0^(9) ~ |g 2 lng| 1//2 , agrees with RPA 
calculations however, the effect of g±, which is a short-range exchange contribution, 

are usually neglected in those calculations. 

The spin part of the Hamiltonian does not involve the long-range part of the interaction 
and only depends on the backward scattering amplitude g\. For repulsive interaction, the 
long-wavelength spin excitations then are described by a Hamiltonian similar to the first 
term in (||), giving rise to collective spin oscillations with u) a (q) = u a \q\, and spin wave 
velocity u a = VpJl — g\. Together with the charge oscillations @, these excitations are 
the complete spectrum of the model. 

The bosonization method makes the calculation of correlation functions rather straight- 
forward. Here, the charge-charge correlations are of particular interest. Using the expression 

p(x) = -(V2/n)d x (j) p (x) 

+ J_ e 2ik F x e -iV2M*) cos[V2<f) a {x)} + cst.e 4ikFX e- iV * Mx) + h.c. (5) 
2ira 

the evaluation of the charge correlation function reduces to the calculation of averages of 
the type 

1/2 

(1 — cosga;) ~ c 2 Vln x , (6) 

with C2 = \J {\ + gijTtv f / e 2 ■ One thus obtains 

(p(x)p(O)) = A 1 cos(2k F x) exp(— c 2 vln~x)/x + A 2 cos(4/c^x) exp(— 4c 2 V / hix) + ... , (7) 

where v4 lj2 are interaction dependent constants, and only the most slowly decaying Fourier 
components are exhibited. The most interesting point here is the extremely slow decay 
(much slower than any power law!) of the 4kp component, showing an incipient charge 
density wave at wavevector 4kp (instead of the usual 2kp of the Peierls instability). This 
slow decay should be compared with the case of short-range interactions, where the 2kp 
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dq l + h 

q l-g 1 + 2V(q)_ 



and 4kp components deacay as with the power laws x~ x ~ p and x~ p , respectively, with 
an interaction-dependent constant K p The 4k p oscillation period is exactly the 

average interparticle spacing, i.e. the structure is that expected for a one-dimensional 
Wigner crystal. Of course, because of the one-dimensional nature of the model, there is 
no true long-range order, however, the extremely slow decay of the 4k p oscillation would 
produce strong quasi-Bragg peaks in a scattering experiment. It is worthwhile to point out 
that this 4kp contribution arises even if the Coulomb interaction is extremely weak and 
depends only on the long-range character of the interaction. On the other hand, any 2k p 
scattering is considerably weaker, due to the 1/x prefactor in (^) which has its origin in the 
contribution of spin fluctuations. 

Other correlation functions are easily obtained. For example, the spin-spin correlations 

are 

(S(x)-S(O)) ^ B lC os{2k F x) exp{-c 2 Vl^x)/x + ... (8) 

where there is no 4kp component. On the other hand, correlation functions that involve 
operators changing the total number of particles (e.g. the single particle Green's function) 
decay like exp[— est. (In x)~ 3 / 2 ], i.e. faster than any power law. This in particular means that 
the momentum distribution function and all its derivatives are continuous at kp, and there 
is only an essential singularity at kp. The calculations are also straightforwardly generalized 
to finite frequency and temperature |l|,@], however the rather complicated formulae are not 
of immediate interest here. 

The presence of metallic screening changes the above behaviour: a finite screening length 
£ s would lead to a saturation of V(q) for q — > at 2e 2 m(£ s /<i). One than would have, for 
x > £sj power-law decay of the type discussed above for short-range interactions, with 
K p wl/ \/ln£ s . On the other hand, if the interaction potential decays more slowly than 1/r 
(a rather hypothetical case), the integral (|6|) remains finite for x — > oo, and therefore there 
then is real long-range order of the Wigner crystal type. 

It is instructive to compare the above result ([?]), obtained in the limit of weak Coulomb 
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interactions, with the case of strong repulsion (or, equivalently, heavy particles). The con- 
figuration of minimum potential energy is one of a chain of equidistant particles with lattice 
constant a, and quantum effects are expected to lead only to small oscillations in the dis- 
tances between particles. The Hamiltoniam then is 

# = Ef^ + ^E V"(ma) («, - u l+m f , (9) 

l l^m 

where u\ is the deviation of particle / from its equilibrium position. In the long-wavelength 



limit, the oscillation of this lattice have energy uj(q) = \J2/ {ma)eq\ ln(ga)| 1//2 . The most 
slowly decaying part of the density-density correlation function then is 

47T 



(p(x)p(O)) ~ cos(27rx/a) exp 



rv/hT 



x 



(10) 



(2me 2 a) 1 /2 

Noticing that k F = 7r/(2a), one observes that the results (|7|) and ([UJ) are (for g\ = 0) 
identical as far as the long-distance asymptotics are concerned, including the constants in 
the exponentials. Eq. (|7|) was obtained in the weak interaction limit, whereas ( [T0| ) applies 
for strong Coulomb forces. Similarly, the small-g limit of the charge excitation energies 
is identical. We thus are lead to the rather remarkable conclusion that the long-distance 
behaviour of correlation functions is independent of the strength of the Coulomb repulsion, 
provided the interaction is truly long-ranged. 

In recent experiments, one-dimensional structures with two partially filled subbands have 
been investigated J7|. If only the long-range part of the Coulomb interaction is considered, 
the appropriate generalization of the model (§) to that case is described by the Hamiltonian 

i=0,l 

+— J dxdx'V{x - x')d x {4> pfi + (j) Pt i)d x i((t) pfi + <p pA ) , 

where Vi is the Fermi velocity of band i, and IT^j are the charge fields of band i. In the 
long-wavelength limit the charge oscillation eigenmodes have energies 



LO + (q) = \q\\j2(v + Vi)V(q)/ir , uj-(q) = v^i|?| 



The uj + mode represents in phase oscillations of the two bands and has the typical \q 2 lngl 1 / 2 
behaviour of one-dimensional plasmons, whereas the o;_ mode is an out-of-phase oscillation. 
In addition there are two spin modes, at energies t>o,i|<z|- If the various possible interaction 
processes involving momentum tranfer of order 2kp,i are included, the velocities of these 
modes are renormalized, similar to the effect of gi in the one-band model above. 
For the density correlations of the two-band model I find 

(p(x)p(O)) = C 2 cos[4(/c Fi0 + k F)1 )x] exp(— 4c 2 Vlnx) 
+x~ K ^2 Ci cos(4k Fti x) exp(— q n 



x) 

i=0,l 



where now c 2 = yjn(v + t>i)/e 2 , q = Av 2 c 2 /(v + vt) 2 , K = A^/v^/(v + t>i), and again 
only the most slowly decaying Fourier components are exhibited. The most slowly decaying 
part of (|T0| ) (q = 4(A^ )0 + &f,i)) is again the one that corresponds to Wigner-crystal type 
ordering, e.g. the electrons order approximately equidistantly. In fact, this type of ordering is 
determined only by the uj + mode, whereas all other Fourier componets contain contributions 
from the mode, which lead to power law decay. 

To compare the present results with experiment J7| one can first notice that, provided 
that 2kpd < 1 and including the background dielectric screening, one has g\ < 0.2, and 
consequently to within a few percent u a = vf, i-6. the triplet spin mode ("SDE") is expected 
at vpq, as experimentally observed. Further, in the experimentally accessible range q < 
0.2kp the plasmon energies found here are indistinguishable from RPA results, and thus the 
present results provide a good fit to the experimental plasmon dispersion. 

More difficult to explain is the extra feature which has been interpreted as an electron- 
hole continuum ( "SPE" ) [|7| : in fact in the present model with its linear electron dispersion 
relation, there is no such continuum (and it would not exist in an RPA calculation either). 
However, the model offers an alternate possibility: together with the triplet spin mode, 



there is also a singlet mode JT1] . The existence of the singlet mode is a consequence of spin- 
charge separation in one-dimensional fermion systems, and in particular it is degenerate 
with the triplet mode. This mode can be found e.g. in energy density correlation functions 



(as opposed to the spin mode, which appears as a pole of the spin density correlation 
function), and therefore also is expected to be seen in the polarized Raman spectra. This 
interpretation requires the SDE and SPE features to appear at the same energy, which seems 
to be consistent with the results published in ref. J?|. It is noteworhty that, if correct, this 
interpretation would mean that these results constitute the first direct spectroscopic evidence 
for the existence of individual spin- 1/2 objects ("spinons"): the existence of degenerate 
triplet and singlet mode implies that they are both build up from non-interacting spin- 1/2 
excitations. 

One might argue that the existence of the particle-hole continuum is due to effects of 
band curvature, which is neglected in the present model. RPA calculations including band 



curvature certainly predict both a plasmon and a particle-hole continuum P,[T0||. However, 
within RPA the total spectral weight for the continuum is about two orders of magnitude 
smaller than that of the plasmon, whereas in ref. |7j plasmon and SPE have comparable 
weight. Moreover, in exactly solved one-dimensional models like the Hubbard model || , 
one finds both a plasmon-like collective mode and the singlet mode discussed above, but no 



separate particle-hole continuum |12[|. There thus seems to be little theoretical evidence in 



favor of an interpretation of the SPE feature in terms of a particle-hole continuum. 

The nearly long-range Wigner crystal type order should have important consequences 
for transport properties: in fact, in the presence of disorder, a classical charge density wave 



(which has real long-range order at T = 0) becomes disordered fl3 |, with a "pinning length" 
describing the decay of spatial correlations given by 

U « {{V/v F ) 2 n)-^ , (11) 

where n is the density of impurities, and V the Fourier component of the impurity potential 
at the wavevector of the CDW (4&i? in our case). Inclusion of quantum effects in systems 
with short range interaction only leads to corrections to the exponent 1/3 in ( |TTD ||14|| . 
Following the same arguments, I expect (|TT|) to be valid for the Coulomb system too (up to 
logarithmic corrections), i.e. as far as low-frequency phenomena are concerned, the system 
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of electrons interacting with Coulomb forces behaves like a classical charge density wave. 
In particular, all the unusual dynamical properties associated with nonlinear transport in 
CDW systems should also occur in the one-dimensional electron system. 

At finite temperature, thermal agitation can become sufficiently strong to depin a CDW. 
In the present case, this is expected to happen when the thermal correlation length, in 
the absence of impurities, given by u p (1/£t) ~ T, becomes shorter than £ pin . 

Nonlinear current-voltage relations characteristic of CDW transport have been observed 



in one-dimensional semiconductor structures |15[ . In these experiments, strong variations of 
the linear conductance with carrier density have been interpreted in terms of the Coulomb 
blockade |16| |, implying the existence of a pair of strong impurity potentials. It seems 
certainly conceivable that the nonlinear behaviour should be due to collective motion of a 
pinned CDW, possibly in parts of the sample where there is only a weak random potential. 
The fact that the order of magnitude of the nonlinear conductance is independent of the 
value of the linear conductance suggests that different mechanisms are involved in the two 
phenomena. The present calculation then demonstrates that Coulomb interactions do not 
have to be particularly strong to create nearly classical CDW type behaviour. It would 
clearly be interesting to investigate the possible interplay between Coulomb blockade and 
CDW-like behaviour. 

In conclusion, using the bosonization technique a consistent microscopic picture of the 
excitation spectrum and the correlation functions of a one-dimensional electron gas inter- 
acting with long-range Coulomb forces has been obtained. The density correlations are 
those of a nearly perfect one-dimensional Wigner crystal. The results are in agreement 
with excitation spectra observed in one-dimensional semiconductor structures and provide 
a possible explanation for nonlinear transport properties. 
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